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William L. Cook 
Goddard Space Flight Center 


ABSTRACT 

The dynamic mode shapes and frequencies of a circular ring 
are obtained using finite element techniques. The results obtained 
from models making use of half symmetry and of quarter symmetry 
are compared with those for the complete ring and with the exact 
solution. 

The half-ring model provides the same solution as the complete 
ring. The quarter-ring model, however, must be subjected to two 
separate sets of boundary conditions in order to obtain the complete 
set of mode shapes. 

This same quarter-ring model is then subjected to a third set 
of boundary conditions to obtain the same solution as an infinitely 
long "corrugated" wire. 
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FINITE ELEMENT MODELING OF A CIRCULAR RING 
USING HALF AND QUARTER SYMMETRY 

by 

William L. Cook 


1. INTRODUCTION 

In the finite element modeling of symmetric structures, a saving in com- 
puter time can often be obtained by determining the response of a segment of 
the structure constrained by appropriate boundary conditions. 

A class of such structures is characterized by radial symmetry (e.g., cylin- 
drical shells, spherical shells, circular disks and circular rings.) The present 
investigation demonstrates a method of approach in applying symmetry-oriented 
modeling techniques to structures possessing radial symmetry. 

Computations were performed on the IBM 7094 computer using the Martin 
Company SB038 program for the analysis of finite element models by the force 
method. 1, 2 


2. ANALYTICAL SOLUTION 

The in-plane vibrations of a circular ring are discussed by Timoshemko in 
Vibration Problems in Engineering , p. 425 ff. 3 

The frequency of pure radial vibration is: 



All displacements are in the radial direction, and hence the motion is de- 
scribed by a circle of periodically varying radius. 
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The frequencies for longitudinal (circumferential) vibrations are given by 


P ‘ = \^ r/T7TT 

where i denotes the number of wave lengths to the circumference. All displace 
ments in this case are in the circumferential direction. The modal frequencies 
for flexural vibrations are: 


i 2 (1 - i 2 ) 2 
(1 + i 2 ) 

u. = a i cos 14> = radial displacement 

v £ = Sl l /i sin i<t>= circumferential displacement in which 
a A is a constant coefficient. 

I = area moment of inertia of cross section 

g = acceleration of gravity 

E = Young's modulus 

r = radius of the ring 

y = specific weight of the ring 

A = cross sectional area 

3. PARAMETER VALUES 

The following values of the parameters were chosen for present investiga- 
tions: 

r = 10.0 in. 

E = 1.0 x 10 lb/in. 


p, =,/** ( — 

y \ Ar 3 4 , 


where the mode shapes are described by: 
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g = 386.06 in./sec 2 


A = 0.19634954 in. 2 
y = 25.9382231 lb/in. 3 
I = 0.0654498467 in. 4 

The values for the last three parameters were chosen so as to provide conven- 
ient values for the grid point masses and element compliances in the model. 


4. SYMMETRY METHODS 

The first finite element model to be investigated is a representation of the 
complete ring, and from this model the entire set of radial, flexural and circum- 
ferential mode shapes and associated frequencies are obtained. 

It will be demonstrated in this report that the same information may be ob- 
tained by modeling only a portion of the structure, constrained by appropriate 
boundary conditions, which results in an appreciable savings in computer time. 

The boundary conditions on the segment are chosen so as to be compatible 
with the deflections experienced by the complete structure in each mode shape 
desired. Hence, any modes whose deflections are not compatible with a given 
set of boundary conditions cannot be obtained from that set. It is obvious from 
the statements above that symmetry methods are applicable only in cases where 
the nature of the mode shapes has been predetermined. 

Consider the two representative mode shapes of the ring shown in Figure 1. 
These are the flexural modes with i equal to 3 and 4 respectively. The radial 
displacements, u ; , are indicated by a solid line, and the circumferential displace- 
ments Vj , by a broken line. 

Assuming the half-symmetry model to lie in the region 0 5 4> - 77 , it is necessary 
to prescribe the boundary conditions at 4> ~ 0 and at <p=Tr . Refering to the modes 
shown in Figure 1 it will be noticed that the slope (du/d^) and the circumferential 
displacement (v) are zero at both points. It can be shown that this condition is 
also true for the circumferential modes. Hence, the half-symmetry model con- 
strained such that these boundary conditions are satisfied may be expected to 
produce the entire set of mode shapes. 
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In the case of a quarter-symmetry model in the region boundary 

conditions must be determined at < t > = 0 and at 4 > =tt/2. As seen in Figure 1 the 
third flexural mode satisfies the condition that the slope (du/d i) and the circum- 
ferential displacement (v) are zero at 4>= 0, and the radial displacement (u) is 
zero at <£= tt/2. This set of boundary conditions will produce all the odd-numbered 
flexural and circumferential modes. 

The fourth flexural mode (as well as all the remaining modes) may be ob- 
tained by constraining the slope and the circumferential displacement to be zero 
at both ends. 

Great care should be taken that all components of displacement are investi- 
gated in the determination of the boundary conditions. For instance, from the 
fact that the radial displacement at <t> = rr/2 is identically zero for the odd num- 
bered modes, it might mistakenly be infered that a node exists at that point (i.e., 
that both u and v are zero.) 

It will be shown that the same quarter- symmetry model discussed above, 
when constrained in just such a manner so that du/d(p and v are zero at ^=0, 
and both u and v are zero at 4> - 77/2, will no longer represent a circular ring. 

The mode shapes in this case will be those of a structure formed by an anti- 
symmetric extension of the segment about =77/2. 


5. FINITE ELEMENT MODELS COMPLETE RING 

The model for the complete ring consists of 32 grid points equally spaced 
about a circle of radius 10.0 inches, connected by 32 primary tension elements 
and 64 two-digit bending elements, 61 of which are primary (Figure 2). 

The structure is supported by 3 reaction elements which are later removed 
for free-body analysis. 

The values for the element compliances and grid point masses are shown 
in Table I. 


Half Symmetry 


The half-symmetry model consists of 17 grid points evenly spaced about a 
semi-circular arc of radius 10.0 inches, connected by 16 primary tension ele- 
ments and 32 primary 2-digit bending elements (Figure 3). 
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Table I 

Compliance and Mass Data for Complete Ring 


DIAGONAL COMPLIANCE 


Element Numbers 

Formula 

Value 

1-32 

C = L/AE 

1.0 x 10" 6 

33-96 

C= L/3EI 

1.0 x 10 -6 


COUPLED COMPLIANCE 


Element Numbers 

Formula 

Value 

33 to 34, 35 to 36, 
... 95 to 96 

C = L/6EI 
1 

0.5 x io" 6 


MASS 





















UNDERLINED ELEMENTS ARE REDUNDANT 

INDICATES REMOVAL PRIOR TO DYNAMIC 
ANALYSIS 

O (O BENDING ELEMENT 

— MOMENT REACTION 
AXIAL REACTION 



Figure 3-Model for Half Symmetry 



Axial reactions 49 and 52 prevent circumferential motion at the ends of the 
ring, and moment reactions 51 and 53 prevent rotational motion at these points. 
Axial reaction 50 is necessary only for stability and is removed prior to deter- 
mination of the mode shapes. 

The compliance and mass values for this case are shown in Table II. The 
reactions in all of the models have zero compliances. 


Table II 

Compliance and Mass Data for Half Symmetry 
and for S- Shaped Curve 


DIAGONAL COMPLIANCE 


Element Numbers 

Formula 

Value 

1-16 

C = L/AE 

1.0 x io -6 

17-48 

C = L/3 El 

1.0 x io 


COUPLED COMPLIANCE 


Element Numbers 

Formula 

Value 

17 to 18, 19 to 20, 
... 47 to 48 

C = L/6EI 

0.5 x 10~ 6 


MASS 


Grid Points 

Mass 

1,17 

5.0/ g 

2-16 

10.0/g 


Quarter Symmetry (Even Modes) 


The circular ring is represented in this case by a quarter segment. 

In order to obtain the radial mode and the even numbered flexural and cir- 
cumferential modes, the ends of the segment are constrained such that rotation 
and circumferential translation are zero. 
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UNDERLINED ELEMENTS ARE REDUNDANT 



GRID POINT 
TENSION ELEMENT 

BENDING ELEMENT 
MOMENT REACTION 
AXIAL REACTION 


Figure 4-Model for Quarter Symmetry Even Modes 



The model (Figure 4) consists of 9 grid points and 28 elements including 8 
tension, 16 bending, and 4 reaction elements (3 of which are primary). 

The compliance and mass values for this case are shown in Table HI. 

Table III 

Compliance and Mass Data for Quarter Symmetry 


DIAGONAL COMPLIANCE 


Element Numbers 

Formula 

Value 

1-8 

C = L/AE 

1.0 x io~ f> 

9-24 

C = L/3EI 

1.0 x 10" 6 


COUPLED COMPLIANCE 


Element Numbers 

Formula 

Value 

9 to 10, 11 to 12, 
... 23 to 24 

C = L/ 6 El 

0.5 x 10~ 6 


MASS 

Grid Points 

Mass 

1, 9 

5.0/ g 

2-8 

10.0/g 


Quarter Symmetry (Odd Modes) 

The same quarter segment used to obtain the even modes is now subjected 
to a different set of boundary conditions to obtain the odd numbered modes (Figure 
5). 


One end of the segment is constrained as before, while at the other end an 
axial reaction (element 28) prevents motion in the radial direction. Reaction 
element 27 is necessary only for static stability and is removed prior to deter- 
mination of the mode shapes. Mass and compliance values are the same as in 
Table III. 
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Quarter Symmetry with Boundary Conditions to Simulate S-Shaped Curve 


The model for this structure is identical to that shown in Figure 5, except 
that reaction element 27 is not removed prior to determination of the mode shapes. 
In other words, all translational motion is arrested at grid point 9. 

The mode shapes obtained using this set of boundary conditions are those of 
a structure formed by an antisymmetric extension of the quarter segment, i.e. , 
of the s-shaped curve described in the following section. 

Mass and compliance values are the same as those in Table III. 


S-Shaped Curve 


This model is the same as the half-symmetry model except that the last 
eight grid points are shifted so as to form an S-shaped curve (Figure 6). Mass 
and compliance values are the same as those in Table II. 

The purpose of this exercise is to provide a comparison for the mode shapes 
obtained in the preceding quarter symmetry mode. The extension of that model 
into an antisymmetric geometry results in the curve considered here. 

It should be noted that mode shapes obtained in this analysis also describe 
the motion of the figure obtained by extending the curve infinitely in either direc- 
tion in a symmetrical manner. Such a figure may be described as a "corrugated" 
wire. 


6. COMPARISON OF RESULTS 

The values of the modal frequencies obtained from both the exact and the 
finite element solutions are compared in Table IV. 

The half-symmetry and quarter- symmetry models are seen to yield the same 
results as the full circle model. 

The differences between the finite element results and the exact results for 
the circular ring may be accounted for by the coarseness of the partition in the 
models. 
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GRID POINT 


UNDERLINED ELEMENTS ARE REDUNDANT 

"INDICATES REMOVAL PRIOR TO DYNAMIC 
ANALYSIS 


MOMENT REACTION 
AXIAL REACTION 


O O TENSION ELEMENT 

A 


o 


O BENDING ELEMENT 


Figure 5-Model for Quarter Symmetry Odd Modes 


KEY: 


UNDERLINED ELEMENTS ARE REDUNDANT 

INDICATES REMOVAL PRIOR TO DYNAMIC 
ANALYSIS 


GRID POINT 


TENSION ELEMENT 



BENDING ELEMENT 


MOMENT REACTION 


AXIAL REACTION 


Figure 6-Model for S-Shaped Curve 



Table IV 

Modal Frequencies (Rad/Sec) 
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